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Abstract 

Using group theory of mixing to examine all finite subgroups of SU{3) with an order less than 
512, we found recently that only the group A(150) can give rise to a correct reactor angle ^13 
of neutrino mixing without any free parameter. It predicts sin^ 2^13 = 0.11 and a sub-maximal 
atmospheric angle with sin^ 2^23 = 0.94, in good agreement with experiment. The solar angle 612, 
the CP phase 5, and the neutrino masses mj are left as free parameters. In this article we provide 
more details of this case, discuss possible gain and loss by introducing right-handed symmetries, 
and/or valons to construct dynamical models. A simple model is discussed where the solar angle 
agrees with experiment, and all its mixing parameters can be obtained from the group A(600) 
by symmetry alone. The promotion of A(150) to A(600) is on the one hand analogous to the 
promotion of 53 to 84^ in the presence of tribimaximal mixing, and on the other hand similar to 
the extension from to 5*4 in that case. 
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I. INTRODUCTION 



Before the reactor angle ^13 was successfully measured in the recent past, the neutrino 
mixing data were consistent with a zero ^13 and a maximal atmospheric angle ^23- These 
two are explainable by the leptonic symmetry S^- This symmetry group is generated by 
two unitary matrices, F = diag(l, w, w^), where cu = exp(27ri/3), and G with = 1, 
det(G') = 1, which possesses the eigenvector (0, 1, —1)"^ with eigenvalue +1. This symmetry 
predicts a zero reactor angle and a maximal atomospheric angle, but leaves the solar angle 
9i2 undetermined. 

The solar angle can also be explained by symmetry if we enlarge the group S3 by in- 
corporating another unitary generator G' which commutes with G, satisfying G'^ — 1, 
det(G") = 1, and possessing an eigenvector (1, 1, 1)^ with eigenvalue +1. The order of the 
group Q generated by F, G, G' must then be an even multiple of 6, the order of S^. 

The simplest non-abelian group with order 12 is ^44, but it contains only even permuta- 
tions so it does not contain 5'3 as a subgroup. The simplest non-trivial group that contains 
S2, is 5*4, with order 24. Moreover, it also contains the G' above, so it predicts the correct 
solar angle. The resulting neutrino matrix has the well-known tribimaximal form. 

Now that the reactor angle is no longer zero and there are indications that the atmospheric 
angle may not be maximal, the question is whether we can still find a symmetry to explain 
the mixing data. In a recent work, we searched all the finite subgroups of SU{2>) with an 
order less than 512, and discovered that the only group that can explain the reactor and 
atomospheric angles by symmetry alone is A(150) [1]. It is generated by the same F as 
S2,, but a different G. It also leaves the solar angle free. In order to have a hope to have 
the solar angle also explained by symmetry, we need to enlarge the group by incorporating 
another order-2 generator G' which commutes with G. As before, the extended group Q must 
have an order which is an even multiple of 150, the order of A(150). Again no non-trivial 
order-300 group contains a A (150) subgroup, but the order-600 group A (600) does, and in 
fact that is the only order-600 group that does. The symmetry of A(600) also contains the 
same G' before which gives rise to a solar angle that agrees with experiment. In this way 
the pair (A(150), A(600)) is analogous to the previous pair (6*3, 5*4) before the reactor angle 
was measured. 

It is interesting that A(600) is the smallest of the three groups that yield a full Z2 x Z2 
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residual symmetry in the neutrino sector, among all the finite SU{3>) subgroups of orders 
< 1536 [2]. It also shows up in the method using von Dyck groups [3, 4]. 

In this article we discuss the groups A(150) and A(600) and their predictions [5]. The 
mathematical properties of these groups are summarized in Sec. II. Given a symmetry group, 
there are three ways to apply its symmetries to neutrino physics, each with an additional 
assumption than the previous, but potentially (though not necessarily) also an additional 
gain. We will discuss these three in turn. The simplest is to apply symmetry only to the 
left-handed fermions, to obtain information on the effective left-handed mass matrices. This 
is the method used in the recent survey [1] and shall be referred to as the 'left-handed 
symmetry' method. The next simplest is to apply symmetry to both the left-handed and 
the right-handed fermions, to obtain information also on the Dirac and Majorana mass 
matrices. We shall refer to this method as 'both-handed symmetry'. Without imposing the 
additional right-handed symmetry there is no way to get the Dirac or the Majorana mass 
matrices. The most involved but the most well-known method is to construct dynamical 
models. This requires the introduction of additional scalar fields called valons to account 
for the spontaneous breaking of symmetry. These three methods will be applied to A (150) 
and their relative merits will be discussed in Sees. Ill, IV, V respectively. 

No matter which of these three methods is used, we can account for the correct reactor 
and atmospheric angles as long as A (150) is the symmetry. The other neutrino mixing 
parameters, as well as the fermion masses, are hidden in the adjustable parameters in the 
mass matrices. For left-handed symmetry, there are just enough parameters to fit all the 
experimental quantities. For both-handed symmetry, the additional freedom in assigning 
right-handed symmetry increases the number of parameters, but except for special cases 
such as the c — example in the text, the additional parameters do not yield additional 
information on the mixing parameters nor the masses. Nevertheless, when dynamical models 
are built, most of these parameters in both-handed symmetry will acquire a physical meaning 
in terms of Yukawa coupling constants. The groups A(6n^) can also be incorporated into 
grand- unified-theory models [6]. 

In principle, a fit to the experimental data can determine the parameters in left-handed 
symmetry, but in practice this is impossible because CP-phase and one of the neutrino 
masses are unknown. The best one can do then is to produce simple models, hopefully with 
simple choice of parameters, that can account for all the known experimental quantities. 
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Such a model will be given in the text, first in Sec. Ill, subsequently followed up and refined 
in Sees. IV and V. 

In Sec. VI, the left-handed symmetry will be applied to A(600), and in Sec. VII, a 

summary of the results will be given. There is also an appendix at the end to explain real 
representations for the group A(150). 

II. A(6n2) 

This series of finite subgroups of SU{3) have been studied in various places [7]. In this 
article we will follow the information given in the mathematical software GAP [8] . 

A(6n^) can be described by the semidirect products of a direct product of cyclic groups. 



The generators of the two cychc groups are /4 and respectively; the generator of 
is /2, and the generator of Z2 is /i. The orders of /i, fs, /2, fi are respectively n, n, 3, 2. 

The direct-product of the two Z„'s tells us that /4 and /s commute. The semidirect with 
Z3 implies that /2/3/2~^ and f^fit^^ are monomials of /s and /4. The semidirect product 
with Z2 implies that /i/s/f^ and fififi^ are monomials of /2, /a and /4, while /i/2/i~^ is 
a monomial of /2. 

With these commutation relations, one can rearrange any monomial of the four /j's, 
namely any group element of A(6n^), into the canonical form /4*/|^/|^/f\ with 64, 63 rang- 
ing between and n — 1, 62 between and 2, and ei either or 1. The total number of 
independent elements of the group can easily be counted from the canonical form to be 6n^. 

In this article we are interested in the case n = 5 and n = 10. We will discuss the details 
of these two groups separately below. 

A. A(150) = (Z5 X Z5) >4ZsXsZ2 

The nontrivial commutation relations of the generators are 



A(6n') = {Zn X Zn) XZ3X Z2. 



(1) 



72/4/2 ^ — hfs, /2/3/2 ^ — /I/I) 



/1/4/1 



-1 
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This group has 13 conjugacy classes, Ci to C13, and 13 irreducible representations, IRl 
to IR13. Its character table is given in Table 1, with various notation explained in Table 
2 and below. The GAP notation /A in Table 1 means the complex conjugate of A, and 
similarly /B,/F, and /G. The notation *C is obtained from C by changing \/5 to — \/5, 
and similarly *D and *E. 

The designations of the IRs in Table 1 are those in the IrreducibleRepresentations com- 
mand of GAP. They differ from those used in the command CharacterTable, which are 
denoted as X.n in GAP, with n varying from 1 to 13. 

The generators /i,/2,/3,/4 belong to classes (^2,^3,^4,^5 respectively. The only class 
with order 2 is C2. No two elements within this class mutually commute. 

As usual, the dimension of an IR is its character for the identity class Ci. The charac- 
ters X of the 3-dimensional representations IR4 to IRll are complex, with the pairs (4,10), 
(5,11), (6,8), (7,9) being complex-conjugates of each other. The other representations, IRl 
to IR3 and IR12, IR13, have real characters, so they are either real or quaternionic rep- 
resentations. Their Probenius-Schur indices X]ggA(i50) x(fl'^)/150 are all +1, so all of them 
are real representations. The representations given by GAP for IRl and IR2 are explicitly 
real, but not those for 1R3, IR12 and IR13. In Appendix A we show how to convert the 
GAP-representations into representations that are explicitly real. 
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Table 1. Character Table of A(150) 
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Table 2. Abbreviations used in Table 1, where 77 = e^'^*/^ 



Using Table 1, we can explicitly work out the Clebsch-Gordan series. The result is given 
in Table 3. Since IRi x IRj is the same as IRj x IRi, only half of the table is explicitly 
shown. For example, the products of IR12 and IR12 is given by (with the prefix IR omitted) 

12 X 12 = 1 + 2 + 3 + 3 + 6 + 7 + 8 + 9 + 12 + 12 + 13. (3) 

Prom Table 2, it can be checked that there are 36 states both on the lefthand side and the 
righthand side of (3). 
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Table 3. Clebsch-Gordan series of A (150) 



B. A(600) = {Zio X Zio) XI Z3X Z2 

Since Ziq — Z^x Z2 (a prime is used to distinguish this Z2 from the one in the semidirect 
product), its generator is related to the Z5 generator C4 and the Z2 generator by 

C4-/4^ di = fi\ U = Ci^di. (4) 
Similarly, for the other Ziq in the direct product, 

C3 = /3', ds^fs^ fs^Cs'ds. (5) 
The nontrivial commutation relations between fi, f2, C3, C4 are the same as in (2), namely, 

1204/2^ = C4C3, /2C3/2"^ = clcl, 

fiC4fr' = 44, hcsfr' = cs, fMr' = fl (6) 
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Together they generate the subgroup A(150) of A(600). 
The other nontrivial commutation relations of A(600) are 

fidif^^ = d^da, /i4/r^ = ds. (7) 

The group formed by {/i,/2,c?3} (or {/i, 72,^4}) is 5*4, thus A(600) contains S4 as a sub- 
group. The subgroup formed by {f2,d3} (or {/2,c?4}) is A4, and the subgroup formed by 
{fijds} (or {/i,c?4}) is the Klein four-group Z2 x Z2. [17] 

Substituting (4) and (5) into (8), we get the nontrivial commutation relations with the 
Zio X ZiQ generators to be 

72/4/2 ^ = /l/s) 12/3/2 ^ = flfi^ 

/1/4/r^ = /4/3 ' /1/3/r^ = fa- (8) 

The group A(600) has 33 classes and 33 IRs, as listed in Tables 4 and 5. We have no 
need for the character table nor the Clebsch-Gordan series in this article, so they will not 
be given here. 

The classes containing order-2 elements are Cg a-nd C24- Unlike A(150), it is possible to 
find two mutually commuting order-2 elements in this group. The significance of this remark 
will be explained in Sec. VI. 
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Table 4. Conjugacy classes of A(600) 



8 



IR 
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11233333363333633 
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3363333666666666 



Table 5. Irreducible representations of A(600) 



III. LEFT-HANDED SYMMETRY FOR A(150) 

In this section we first review the general theory of mixing based on a symmetry of 
effective left-handed mass matrices (to be abbreviated as 'left-handed symmetry') [9]. We 
will also provide much more detailed information on the specific case of A(150). 

Every neutrino mixing matrix U has a natural Z2 x Z2 symmetry generated by 

Gi = ±Uiul =F Mjw] ± Ukul, (i = 1, 2, 3) (9) 

where Ui,U2, are the three columns of U, referred to as mixing vectors below, and {i,j, k) 
is a permutation of (1,2,3). It is easy to see that is unitary, = 1, and GiGj = iG^. 
Gi is a symmetry of the Major ana neutrino mass matrix M^, — in the basis where the 
left-handed effective charged-lepton matrix Mg := M^Me is diagonal. This is so because 

GjM^Gi = (10) 

follows from Mj, = ^iTriiUiuf , with rrii being the masses of the active neutrinos with 
Majorana phases. 

This symmetry in the neutrino sector is accompanied by a symmetry in the charged- 
lepton sector. Since the Mg is diagonal, every unitary diagonal 3x3 matrix F satisfying 
F^Mf.F — Me is a symmetry. If 1, then F generates the said Z„ symmetry. F cannot 

be identical to any Gj, otherwise the mixing column Ui would have two zeros, contrary to 
experiments. 

and Mg are the (effective) left-handed mass matrices mentioned earlier. 
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So far everything is general. The group theory of mixing is based on the assumption 
that F and at least one Gi (hereafter referred to simply as G) originate from an unbroken 
nonabelian symmetry Q, and remain a symmetry after the spontaneous symmetry breaking, 
although the other members of Q may not. From now on we shall assume ^ to be a finite 
group. 

The minimal symmetry group Q consistent with this assumption is generated by F and 
G. Conversely, given any we can construct mixing vectors v (one of the Ui above) in the 
following way. 

Collect from Q all possible pair of elements (F, G) belonging to all possible 3-dimensional 
irreducible representations in which F and G are unitary. In order to ensure Mg to be 
diagonal when F is, we require the three eigenvalues of F to be distinct, hence its order 
n must be at least 3. G should be of order 2 because G^ — 1. Let v — (f i, ^2, f^s)^ be 
the normalized eigenvector of G with eigenvalue — Tr(G) (= ±1) in the F-diagonal basis, 
then i> is a possible mixing vector and it is uniquely determined from G. In particular, if 
it is the mixing vector for the third column of the neutrino mixing matrix [/, then \vi\ — 
I sin ^13 1, 1^2 1 = I cos sin ^23 1 , and \v^\ = | cos ^13 cos ^23 1 • Since the CP and Majorana 
phases are unknown, many distinct v's can give rise to the same \vi\ which is all that can 
be measured. Moreover, by rearranging the entries of the diagonal F, the components of Vi 
are similarly rearranged. We shall refer to these distinct i^'s that give rise to the same \vi\ 
with some rearrangement as equivalent mixing vectors. For experimental comparisons they 
need not be distinguished. 

In addition to the G whose mixing vector v fits the third column of U, if we can find 
another G' commuting with G whose mixing vector v' fits the first or the second column, 
then the unitarity of U ensures the remaining column to agree with experiment. In this way 
we could get a purely symmetric explanation of all the mixing parameters, without invoking 
any free parameter. In other words, we have then a full Z2 x Z2 rather than a single Z2 
symmetry in the neutrino sector. A necessary condition for this full symmetry to happen is 
to have two mutually commuting order-2 elements G, G' in the group. As remarked in the 
last section, this can never happen in A(150), but it does happen in A(600). 

There are clearly many (F, G) pairs for each group. The choice in G is somewhat limited 
because it must belong to a conjugacy class of order 2, but F can come from essentially 
everything else except the identity. Fortunately, it is sufficient to take only one F per 
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conjugacy class of order > 3, provided we choose G to be all possible elements of order 2. 
This is so because the pair (F, G) and the pair {gFg~^, gGg'^), for any g E give equivalent 
mixing vectors v, so other F's in the same conjugacy cannot yield an inequivalent mixing 
vector. 

Let us now specialize to ^ = A (150). As shown in Tables 1, it has 13 classes Cj and 
13 irreducible representations (IR). G must be taken from C2, with IC2I = 15, but F can 
be taken from C3 to C13, hence there are 11 possible F's. There are eight 3-dimensional 
irreducible representations (IR4 to IRll), so altogether there are 15 x 11 x 8 = 1320 distinct 
(F, G) pairs. Most of these pairs do not yield an acceptable mixing vector for the third 
column, but many pairs do, all giving equivalent mixing vectors |f| = (.170, .607, .777)"^, 
corresponding to sin^ 2^^i3 = 0.11 and sin^ 2^^23 = 0.94. These pairs occur in IRn, with F in 
Ci, and G being the kth element of €2- Table 6 consists all these possible triplets {n, i, k). 
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Table 6. Origin of equivalent mixing vectors 
We see from Table 6 that the successful F comes from C3, whose GAP-basis expression 
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turns out to be 

/O 1\ 

F= 10 (11) 
\0 1 0/ 

for all the IRs in Table 6. It can be diagonalized by V, yielding V^FV — diag(a;^,a;, 1), 
with (jj := exp(27ri/3), where 



V = 



1 

71 



V 1 



2 UJ 1 

1 1/ 



(12) 



If we permute the columns of V , then the entries of the diagonal form of F are similarly 
permuted. 

The mixing vector v in the GAP basis is the eigenvector of G in that basis, with eigenvalue 
— Tr(G'). The mixing vector v in the F-diagonal basis which yields a possible column of U 
is related to it by f = V^v. 

It turns out that there are only six different G and v pairs, to be labelled by a, 6, c, d, e, / 
below, a comes from the first two rows of Table 6, h the next rows, so on, and / comes from 
the last two rows. They are 
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where 77 := exp(27ri/5). There are four big columns in Table 6. The G's coming from 
columns 1 and 3 have det(G) = —1, and those coming from columns 2 and 4 have det(G) = 
+1. 

The corresponding unnormalized v are 
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The corresponding normahzed mixing vectors v — V^v are 
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Since complex conjugation interchanges uj with u;^, and r^^ with r^^, and since \v\ — \v*\^ we 

conclude that \va\ — — and \vb\ — — \vj\. Moreover, \v}\ is obtained from by 

interchanging the last two rows. Hence all six cases give rise to the same equivalent mixing 

vector whose normalized magnitude is \v\ = (.170, .607, .777)'^. 

The mass matrix Mg = m\ can be obtained from the symmetry condition F^M^F to be 

(a /3 /?*\ 



(16) 



where a is a real parameter and (3 a complex parameter. The masses ml, and can be 

obtained from a, /3r :— ^{/3) and /3i :— to be o; — + V^/3i, a — /3r — y/3/3i, a + I^r. 

Similarly, the neutrino mass = can be obtained from the symmetry condition 

M^G = Mij. Corresponding to the six G^, we get respectively 

/ah c\ / a b c \ /ah 



h ar] cri^ , Mj,6 

\c crj^ f / 
fa b bri^\ 



\bri^ c 



1 4- 2 

ar] cr] 

\c cri^ f J 
[a h 



h f hrf^ 
\ c hrf arj ) 



h f hrf^ 

\ c brj^ ar]'^ ) 
(a b 6r?^\ 



> M.f = 



b frj c 
V6r?2 c / / 



(17) 



where a,b,c, f are complex parameters which can be used to fit the solar angle, the CP 
phase, the three neutrino masses and the three Majorana phases (one of them is an overall 
phase) . 

A. A Simple Model 

Prom now on we will study case a in more detail, and will write Mi,a simply as Mi,. As 
has already been remarked, the four complex parameters a, b, c, / are just enough to fit the 
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neutrino masses and Majorana phases, the solar angle and the CP phase. In this subsection, 
we consider a simple model with c — 0, which turns out to give the right solar angle as well. 
The remaining parameters a, b, f can be used to determine the neutrino masses mj and the 
Majorana phases. 

With c = 0, the mass matrix becomes 

/a b 0\ 

b ar] . (18) 
\0 // 



(M,)o 

It can be diagonalized by the matrix 



/ ^3 Q \ 



V2 



1 
V V2 



1 

/ 



(19) 



resulting in 



/ ar] — br]^ 
l^(f (M,)oWo = 
V 

The neutrino mixing matrix in this case is 





/ 
ar] + br]^ ) 



(20) 



1 



V6 



\ l + r/3 



(21) 



where V is the matrix in (12) used to diagonalize F. 

The third column of C/q is just of (15), which gives the desired reactor and atmospheric 
angles. The second column gives trimaximal mixing, yielding a solar angle ^12 satisfying 



sin^ 2^12 = 0.90, 



(22) 



consistent with the PDG value [10] sin2 2^i2 = 0.95 ± 0.10 ± 0.01. The first column is 
determined from the other two by unitarity of Uq. 

The CP phase can also be read off from the (21), (22), (31), or (32) element of |(C/o)y|. 
The result \s 5cp — 

The neutrino masses are the absolute values of the diagonal elements in (20): 



mi = |a?7 - &77^|, m2 = |/|, 777.3 = 1077 + 677^ 



(23) 
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There are enough parameters to accommodate the measured mass gaps, either with a 
normal hierarchy or an inverted hierarchy. In fact, it can fit any Major ana phase as well. 
For example, if it should turn out that the Majorana phases are all zero, then we would take 
at] =p bi]^ and / to be positive. In a normal hierarchy with mi = 0, then = is the solar 
gap and = 4|6p is the atmospheric gap. We can do so similarly for an inverted hierarchy. 

The charged-lepton mass squares are given below (16) to be ml — a — /3r + y/3/3i, = 
01 — Pr — V^Pi, and ml — a + 2/3ji. Thus — ml — —2\^/3i tells us that /3i is negative, 
a — {ml + + m^)/3 is positive, and 2f3R — m^. — a is also positive. 

IV. BOTH-HANDED SYMMETRY FOR A(150) 

The left-handed mass matrices Mg = MjMg and = MJM^^M^ given in (16) and 
(17) are the most general that yield the correct reactor and atmospheric angles from A(150). 
They have just enough free parameters to fit the unknown quantities: charged-lepton and 
neutrino masses, the solar angle, as well as the CP and the Majorana phases. Left-handed 
symmetry alone cannot determine these parameters, nor can they tell us what the Dirac 
mass matrices Me, Mi, and the Majorana mass matrix Mn are. For the latter we need to 
impose additional assumptions regarding how the right-handed leptons transform. 

Suppose the right-handed charged leptons transform according to a (irreducible or re- 
ducible) representation C, and the right-handed Majorana neutrinos N according to a 
representation C. Let Fc , Gc be the unitary matrices for the residual symmetries in these 
representations. With these symmetries, the Dirac and Majorana mass matrices should obey 
the constraint 

Me = F^MeF, = G^M^G, = G^MnGc, (24) 

where F is given in (11) and G is one of the six equivalent forms in (17). In what follows 
we shall confine ourselves to case a, namely. 





/o 





1\ 






^3 






/ 


-r?3\ 


F = 


1 








, G = - 


^2 












1 




u 


1 


o) 




lo 









V 


/ 
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and 



Me = 



/a 13 f5*\ 

f3* a P 
\I3 a) 



/a b c \ 
My = b arj erf 
\c erf f J 



(26) 



These matrices are taken from IRS, with F coming from class C3 and G from 62- More 
specifically, in terms of the generators fi discussed in Sec. II, 



In order for the left-handed matrices = M^Mf. and Mi, = Mj M]^^ to be given 
by (26), Gc must be real as well as unitary, a requirement which can also be seen from the 
Majorana nature of N. If N is described by a Majorana field ■0, then t/;^ — Ct/j — ■0, where 
C is the charge conjugation operator. To preserve the Majorana nature under an internal 
symmetry transformation ^0 — >■ Rip, it is necessary to have R — R* because ipc — >■ -R*'0c- As 
a result, G — G^ because G is unitary and of order 2, and M^^ transforms like M^: 



With the introduction of right-handed symmetry, we expected more parameters to appear 
in Me, Ml,, and Mjv than in Me and M,,. Nevertheless, when we use the former three to 
calculate the latter two, we must still get back to (26) and the parametrizations there. On 
the one hand, since there are more parameters in the former than in the latter, we expect the 
resulting parameters in the latter to remain independent. On the other hand, the composite 
nature of the latter may produce occasional surprises arising from special composite features. 
A case in point which we will discuss later occurs in the model in Sec. IIIA, where c = 0. 
There are four different right-handed assignments giving the same Mj,, and in three of the 
four cases, 777.3 is forced to be zero as well. 

If future experiment should reveal that 7773 7^ 0, then this is a strong support for the 
symmetry assignment of N to be the fourth. This conclusion is what we gain by imposing 
the right-handed symmetry. 

We will assume fermions of different generations to be distinguished by different family 
quantum numbers, hence if C and C are reducible, the irreducible representations they 
contain must not repeat themselves. Also, included in the latitude of choice of C is the 
number p of A'^'s, which is unknown. We reject p = 1 to retain the possibility of leptogenesis, 
and will confine ourselves here to p = 2 and 3. 



F = /2, G^hflfl 



(27) 



MZ^ = G'^MZ}G. 



(28) 
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A. Charged-Lepton Sector 



There are three right-handed charged leptons, so the dimension of C must be 3. If it is 
irreducible, from Table 1 we see that it must belong to one of IR4 to IRll. If it is reducible, 
it could either be IR3+IR1 or IR3+IR2. 

It turns out that Fc = F ii C is one of IR4 to IRll. Even so, the solution of (24) for 
Me is not the same as solution Me in (26) because Me has to be hermitian but Mg may not 
be. The general solution of (24) for Me turns out to be 

/C ^ r\ 



M,= 



(29) 



T C a 

If it is hermitian, then t = a* and we get back to the form of Me in (26). In any case, 
we can compute Mg = M]Me from (29) and obtain (26) with a = |CP + kP + and 
/3^(*a + a*r + t*C. 

Both Me and Me give the same physics, so how come the former is described by three 
real parameters, a, Pr, Pi, but the latter is expressed in terms of three complex parameters 
C, o", r? This puzzle can be solved by computing V'^M,y to diagonalize Me. One finds its 
three eigenvalues to be complex. Their magnitudes arc the masses and their phases are 
unphysical and can be absorbed into the right-handed charged-lepton fields. 

Next, suppose C is IR3-I-IR1 or IR3-I-IR2. Then using the real representations given in 
(A8), we obtain 

/-I -V3 0\ 
Fc> = \ \/3 -1 
\ 2/ 
for both cases. The symmetry constraint (24) then yield 



(30) 



Me = 



-(2(7 + 0/^3 ((7 + 20/^/3 ((7-C)\/3 



V 



r 



If we compute Mg = MjMe from (31), we get back to (26) with a = 2{\C,\'^ + \a\ 
a\f/2> + |t|2 and /3 = 2(-|C + - a*C)/3 + jrp. 



(31) 

IC + 
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B. Neutrino Sector 

As mentioned before, the representation of a Majorana neutrino must be real. Since 
IR5 is complex, it is not possible to realize M,^ in a type-II seesaw mechanism. There is 
no problem to realize it with a type-I seesaw, as implicitly assumed before, as long as the 
representation C of A?^ is real. The only real representations of A (150) (see Table 1 and 
Appendix A) are IRl, IR2, IRS, IR12 and IR13. The dimensions of IR12 and IR13 are six, 
so we can ignore them if p < 6. This leaves only two possibilities each for p = 2 and 3. 
For p = 2, C =IR3 or IR1+IR2. For p = 3, C is either IR3+IR1 or IR3+IR2. In fact, if 
p < 6, the only other possibility is p = 4 with C=IR3+IR2+IR1, but that contains many 
parameters and sheds no light on the existing physics so we will consider it further. 

We consider these four cases separately below. 



1. p = 2, C= IRS 



From (27), Table 1 and (A8), we get 




(32) 



The solution of (24) is 




(33) 



(34) 



Computing from these — M'^ Mj/M^,, we obtain (26) with 



a 



c 



z{x + yrf){P 



2xyrfQ, b = 2xyP - {x^ + y^r])rfQ, 
Q), f^2z\P-Q), ^ 



a + brj^ 



{x + yri^iP 



Q)- 



(35) 



2. p = 2, C= IR1+IR2 



From Table 1, we get 




(36) 



18 



The solution of (24) is 

/P 0\ fx -XTf 0\ 

M-i= , M.:= . (37) 

Vo Q/ \y yrf z) 

(38) 

Computing from these M^, = M]^ M^^M^, we obtain (26) with 

a = x^P + y^Q, b = rf {-x^P + y'^Q) , 
c = yzQ, f = z^Q, 
a + hri^ = 2y2g. (39) 



3. p = 3, C =IR3+IR1 
In this case, 







1 








Q 


R\ 




( 


X 


y 


' 1 




Gc = 


1 










Q 


p 


R 






-yr}^ 


—xrj^ 


—z 


(40) 




u 





ij 




\R 


R 


s) 




\ 


w 


—wrf 


) 





(41) 

Note that the 3x3 Dirac mass matrix here is a composite, with its first two rows taken from 
the My in (33), and the last row taken from the first row (with a change of symbol) of the My in 
(37). 

Computing from these My = MjM^^My, we obtain (26) with 

a = {x^ + y^rf)P-2xyrf'Q + 2{x-yr]^)wR + w'^S, 
b = 2xyP - {x^r] + y^)ri^Q + 2{-xri^ + y)wR-'u?rf'S, 
c = z{x + yri^){P-Q), f = 2z^{P-Q), ^ 
a + bri'^ = {x + yrl^f{P-Q). (42) 
If we set w = = 5 = 0, we return to (35) as it should. 



I p = 2, C=IR3+IR2 
In this case. 







1 0\ 




fp 


Q 


R \ 
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X 


y 






Gc = 


1 







Q 


P 


-R 


, My := 




-yrf 


—xrj^ 


—z 


■ (43) 
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-ij 




\R 
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(44) 

The 3x3 Dirac mass matrix Mi, here is also a composite, with its first two rows taken from the 
Mjy in (33), and the last row taken from the second row (with a change of symbol) of the M^, in 
(37). 

Computing from these Mjy = MJM^^M^, we obtain (26) with 

a = {x^ + y^rf)P -2xyr}^Q + 2{x + yri^)wR + w'^S, 

h = 2xyP - {x'^r] + y^)ri^Q + 2{xrf + y)wR + w^rfS, 

c = z{x + yrf'){P — Q) + (2zw + u{x + yrf)^ R + uwS, 

f = 2z'^{P-Q)+4zuR + u^S, 

a + hq^ = {x + yrl^)[{x + yrj^){P -Q)+4:wR + 2v?S. (45) 

If we set = = i? = S = 0, we return to (35) as it should. 

C. General Remarks 

In all the cases considered above, both in the chargcd-lepton and in the neutrino sectors, there 
are more parameters than necessary to fix Mg and My in (26), from which all the low-energy 
leptonic physical quantities can be determined. This is specially so in case 4 above. As a result, 
the situation is very complicated and there are many degenerate parameters that cannot be resolved 
by the measured quantities. In the rest of this section, we will concentrate on studying the model 
c = in more detail, where things become a bit simpler and more transparent. 

D. The Model c = 

Consider M^v in (26) with c = 0, a case already considered in Sec. IIIA. Recall that all the 
experimental mixing angles can be obtained with this choice, together with 5cp = 0. We would 
like to know what happens if we also impose a right-handed symmetry in the model. For that 
purpose, let us concentrate on equations (35), (39), (42), (45), and (20). 

It can be seen from these equations that / = often follows from c = 0. This is undesirable 
because phenomenologically m2 = |/| can never be zero, although either mi or may be. In the 
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first three cases, the only way to render c = and / 7^ is to set 

x + OT^ ^ in (35), y = 0in(39), x + yr)^ = in {42). (46) 

In all three cases, this automatically implies ms = |a + 677^ | = 0. This unexpected prediction of 
777-3 = 0) S'lid hence an inverted hierarchy, stems from the right-handed symmetries imposed on 
these three cases, so at least in these cases, the imposition of right-handed symmetry does produce 
additional predictions. 

Beyond that, there are enough free parameters left in each case to fit any 7711 and m2- Specifi- 
cally, for case 1 in (35), 

mi = \a-bri'^\ = \4x'^{P + Q)\, m2 = \ f\ = \2z'^{P - Q)\. (47) 

For case 2 in (39), 

7771 = |a - bri^l = \X^P\, 7772 = I/I = \z^Q\- (48) 

For case 3 in (42), 

mi = \a-br]'^\ = \Ax'^{P + Q) + 8xwR + 2w'^S\, m2 = \ f\ = \2z'^{P - Q)\. (49) 

The last case is more complicated because it involves by far the most number of parameters. 
There seems to be many ways to impose c = and / 7^ 0, and it seems quite possible to have 
7773 7^ with c = 0. For now let us just concentrate on a simple choice of parameters to render 
c = 0, by letting x + yrj'^ = and w = 0. With this choice, 7773 = as well, and as before there are 
enough parameters left to fix 

mi = \a-bri^\ = 2\a\=4\x^{P + Q)\, m2 = \f\ = \2z'^{P - Q) + 4zuR + u^S\ (50) 

in many ways. 

V. DYNAMICAL MODELS FOR A(150) 

It was assumed in the last two sections that the family symmetry Q = A (150) broke down to the 
residual symmetries generated by F and G, leaving it unspecified how to achieve that dynamically. 
In this section we make a further assumption that this is caused by the vacuum expectation values 
of scalar fields coupled to fermions in ^-invariant Yukawa interactions. We will assume these 
scalar fields cf), called valons, to carry only ^-quantum numbers but no Standard Model quantum 
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numbers. The Standard-Model quantum numbers are carried by the usual Higgs field H, so the 
compound field H<p appearing together in the Yukawa terms carries both Standard Model and 
family quantum numbers. Since we are only interested in the family structure here, we will omit 
all spacetime and Standard- Model details, and replace the Higgs field by its expectation value {v). 
In the presence of the compound scalar field, the Yukawa terms have dimension 5, so the coupling 
constant is inversely proportional to some heavy scale A. The Yukawa terms can then be written 
in the form 

E ^5 ^Xa^b {Cc\Bb, Aa) + h.c., (51) 

B,a,b,c 

in which the factor {v)/A has been absorbed into the Yukawa coupHng constant hs- The fermion 
fields and x are assumed to transform according to representations C and A, and the valon field 
according to representation B. The indices o, 6, c are the components of A, B, C, and {Cc\Bb, Aa) 
are the Clebsch-Gordan coefficients needed to render (51) ^-invariant. 

Let gA,9B,9c be an element of Q in representations A,B,C respectively. A basis is chosen so 
that these matrices are unitary, and real if the representation is real (see Appendix A) . That is also 
the basis in which the Clcbsch-Gordan coefficients {Cc\Bb, Aa) are computed. From the covariance 
relation of Clebsch-Gordan coefficients [11] 

Yl i9c)cc'i9B)b'bi9A)a'a{Cc'\Bb',Aa') = {Cc\Bb,Aa), (52) 

a',6',c' 

we see that the Yukawa term (51) is invariant under any ^-transformation 

il^^9Cip, 9AX,4> ^ 9b4>- (53) 

To break the symmetry from Q down to the residual symmetry, we impose vacuum expecta- 
tion values (0^), determined by a ^-invariant valon potential V(0) via the equation of motion 
dV/d(f)^* = 0, where D* is the complex conjugated representation of D. In order to preserve the 
residual symmetry, the invariance (53) must remain true for g = F in the charged-lepton sector, 
and g = G in the neutrino sector. This requires 

FB{ct>^)e = {ct^^)e, Gb{^^). = {^''). (54) 

for every B, where the subscripts e and v indicate in which sector the vacuum expectation values 
apply. Incidentally, since the normalization of (0^) cannot be determined from (54), it is conven- 
tional to take out the energy unit with an arbitrary numerical constant and absorb them into the 
Yukawa constant hs- What remains is usually called the vaccum alignment. 
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In other words, the ahgnment in the charged-lepton sector must be an invariant eigenvector of 

F, and the ahgnment in the neutrino sector must be an invariant eigenvector of G. By invariant 
eigenvector, we mean an eigenvector with eigenvalue +1. 

The mass matrix after symmetry breaking can be read off from (51) to be 

(Me)c„ = Y.^B{<t>f)e{C'c\Bh,Aa), {M,)ca = Y,hB{<t>f)u{Cc\Bh,Aa), 

BM B,b 

{MnU = ^hB{4>b)ACc\Bb,Ca), (55) 

B,b 

where the range of B,b and the value of hs may be different for the three cases. The Q- 
representation of the left-handed isodoublets is A, that of the right-handed charged leptons is 
C, and that of the right-handed neutrinos N is C, with C being a real representation. It follows 
from (52) and (54) that the constraint on the mass matrices so obtained, 

Me = Fl,MeFA, M, = GIm,Ga, Mn = GIMnGc, (56) 

is the same as (24), except that in (24) Fa was abbreviated as F and Ga was abbreviated as 

G. In other words, dynamical models are consistent with both-handed symmetries and may be 
considered as a refinement of the latter. 

Now that we have decided on the desired alignment to preserve residual symmetries, the question 
is whether a potential V, invariant under Q, can be devised to provide such alignments. It turns out 
that every invariant potential does provide alignment solutions which are invariant eigenvectors of 
F, or G, or any other element g of Q. The reason for this will be explained below. 

Let us assume V = V^^^^ + V^^^ to be a polynomial consisting of two parts: V^^^^ of degree 3 
and 4, and V^^^ of degree 2. The quadratic part can be written as 

V('^=-J2^^lc|>^'c|>^, (57) 

B 

where (f)^ = {(j)^)* is (p in the representation B*. Accordingly, the equation of motion can be 
written as 

The unspecified index B in the argument of simply means that this may be a function of 
for several -B's. To obtain {(j)^) in the lowest order, we may simply consider the fields (f)^ in (58) 
to be classical and identify them with (0^). 
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Let gB he an element of Q in representation B. We will now show that if (f)^ is a normalized 
invariant eigenvector of gs, with gB{4>^) = {4>^), then 0-^ = K(f>^ is a solution of (58) for every D 
provided (j)^ is the unique invariant eigenvector of go- The normalization factor k is a constant to 
be determined from (58) in a way to be explained. 

To show this assertion, first note that 

gDQDi<l>'') = QDigBcp"") = QDi<p''), (59) 

so Qi3 is an invariant eigenvector of go, whatever Qd and the normalization of (f)^ are. Suppose 
Qd consists of a number of monomials q]^ of degree 2 and a number of monomials cj^ of degree 3, 
so that Qd = Yli^D + Y^j^D- Since every q\){(l)^) and every c^(0^) is an invariant eigenvector 
of gs, and since this eigenvector is unique, we must have q^i^'^) = ai4>^ and Cj{4>^) = Pj^^ 
for some constants aj and /3j that are determined by the structure of and o'. Since q^ is 
quadratic and cubic, it follows that q'j-){K(f)^) = aiK^4>^ and Cj{4>^) = Pjk^^^, hence Qd{i^^^) = 
Q^i'^^+Sj PiK-^)4>^- Consequently, if k 7^ is chosen to be the solution of the quadratic equation 

= '^aiK + ^PiK'^, (60) 

i j 

then (f)^ = k4>^ is a solution of (58). | 

Let us now denote the that couples to fermions in the charged-lepton sector of (51) by (j)e, 
and the one that couples in the neutrino sector by (py. What we want then is a solution with 
{4'e) ~ {4'^)e a-iid {4>^) = {4>^)y. This can be accomplished if we choose the valon potential to be 
V = Ve(0e) + V;/(0;/), with Ve and Vv being any two ^-invariant potentials. However, if V includes 
an interacting potential Vev{<i>e-, ^v), then like spin-spin interaction, the resulting alignments of {<i)e) 
and {(f)^) would be shifted. To be consistent with everything up to now, we must assume either 
additional 'shaping symmetry' can be found [12] to effectively forbid Vei/, or that Nature provides 
us with a small Vejy so that everything said up to now would be approximately true. 

The dynamical model reviewed above is the most popular method to implement family sym- 
metry, but it requires the presence of valon fields which may or may not exist in Nature. Since 
the Yukawa terms are of dimension 5, the theory is at best a lower-energy effective theory used 
to explain leptonic masses and mixing. The heaviest leptonic mass is mr < 2 GeV, so one might 
expect the valon-degree of freedom to show up around there, but certainly there is no sign of any 
of them up to the present. Moreover, if there are really new degrees of freedom present, one might 
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expect the Higgs coupling to fermion pairs to deviate from the Standard Model [13] , but hitherto 
there is no sign of that either in the 125-GeV Higgs candidate. 

In contrast, either the left-handed or both-handed symmetry approach to the mass matrix does 
not require the presence of valons, though they do leave open the question of how the family 
symmetry is broken to the desired residual symmetry. In this sense they are similar to the texture- 
zero approach to mass matrices, which for example yields the celebrated relation of Cabibbo angle 
in terms of ratio of quark masses [14] , but the origin and the location of the zeros generally remain 
somewhat obscured. 

With the assumption of the existence of valons, the dynamical model does allow a better control 
of symmetry, because it provides a physical explanation for the parameters appearing in M^jM^,, 
and Mn as Yukawa couplings, which in principle can be directly measured if valons do exist. 
Moreover, one might imagine the possibility of turning off certain Yukawa couplings on the grounds 
that the corresponding valons do not exist, thereby obtaining relations between the parameters of 
mass matrices, which may provide relations between physical qTiantitics. 

We do not expect this last possibility to materialize in the charged-lepton sector. As discussed 
in the last two sections, there are three positive parameters present to determine the three charged- 
lepton masses. If a relation between these three parameters exist, it would imply a relation between 
the charged-lepton masses. The only known relation between these masses is the Koidc relation 
[15], but that involves square root of masses that cannot be obtained by this kind of Yukawa 
terms. For that reason we will concentrate on the neutrino sector in what follows. The four cases 
considered below corresponds to the four cases studied in the last section. In these cases, the left- 
handed neutrinos u always transform in representation A =IR5, but the right-handed neutrinos N 
transform according to different representation C in different cases. With p being the number of 
A'^'s, the first two cases have p = 2 and the last two have p = 3. 

Before going into the details of these four cases, it might be useful to provide a brief summary. 
First of all, there are always enough Yukawa coupling constants /ib in every case to provide the 
necessary parameters x, y, z, u, w, P, Q, R, S in the last section. Since some of these parameters are 
redundant in determining the parameters a, 6, c, / of M,y needed to fix the measured quantities, 
some of the Yukawa coupling constants can indeed be set to zero without giving rise to any new 
physical predictions. Next, one might ask whether the model studied in the last two sections with 
c = can be obtained by the absence of some valons, or equivalently, by setting some Yukawa 
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couplings to zero. The answer is no, though c = can be obtained by specific relations between 
two non-zero Yukawa parameters. 

The Clebsch-Gordan coefficients needed to calculate the mass matrix in (55) can be obtained 
with the help of GAP using (52), which states that {Cc\Bc, Aa) is the invariant eigenvector of 
9c ® 9b ® 9a- With that the computation of the Clebsch-Gordan coefficients is reduced to an 
algebraic problem of finding the invariant eigenvectors from the representations of group elements 
given by GAP. 

If {Aa\,{Bh\ are the a, 6 components of two vectors in the representation spaces A and B 
respectively, then 

(Cc| = ^{Cc\Bh, Aa) {Bb\ {Aa\ (61) 

a,b 

is the c component of a vector in the representation space C. In what follows we will list a few 
Clebsch-Gordan coefficients using (61) and the following notation, a is a singlet in IRl, 6 is a 
singlet in IR2, (c, d) is a doublet in IRS, and [e, f, g), {e' , f , g'), {e" , f" , g") are triplets in IRS, 
IRIO, IRll respectively. Subscripts may be used to distinguish two multiplets transforming the 
same way. The product of a multiplet with a always reproduces the multiplet, so we will skip those 
relations. Some other useful ones are 



c 



1 / {l-VS)ee' + il + VS)ff'-2gg' 



dj 2V3 \^_(i + ^)ee' + (-1 + ^3)//' + 255V ' 

'c\ _ 1 /{~l-VS)ee" + {-l + V3)ff" + 2gg"\ 

d) ~ 2V3[(-i + ^)ee" + {-l-VS)ff" + 2gg") 

' c\ I / C\C2 — C\d2 — d\C2 — d\d2 \ 

d J ^ \ — C1C2 — cid2 — diC2 + did2 J 



a = l={ee" + fj" + gg'% 

h = l={ee' + ff' + gg'), 

a = -^{ciC2 + did2), 

b = ■^{Cid2-C2di). (62) 



We will now proceed to study the four cases separately. Since there is really nothing new we 
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can get that way, all that we can do is to relate the parameters x,y, z,u,w, P,Q, R, S of the last 
section to the Yukawa coupling constants hs, and to determine the relations between the /ib's that 
can render c = 0. 



A. p = 2, C=IR3 



Dirac mass matrix 



From Table 1, we know that IR5xIR10 and IR5xIRll both contain IRS, so B in (55) consists 



of IRIO and IRll. The corresponding representation Gb for the residual symmetry is 

/ 0\ 



GiQ — —Gil 



-3 



(63) 



V 1/ 

The alignments obtained from (54) are {(f)^^) = in one case, and an arbitrary linear 

combination of (^^'^") = (77^, 1,0)^ and {(j)^^'^) = (0,0, 1)^ in the other case. The reason why IRll 
has one invariant eigenvector and IRIO has 2 is because det(Gii) = +1 but det(Gio) = There 
are two Yukawa coupling constants hio and hu, but since the combination of the two invariant 
eigenvectors of IRIO is arbitrary, there are effectively three unknown coefficients, which we will 
refer to as /iioai^iofe) S'ld hu. Using (55) and (62), we find that Mj^ is of the form given in (33), 
with 



X=^{hiOax''^'^ + hiix''), y=-l=(/iio„y'0- + /iny'l), z = -^hiob, (64) 



where 



X = r] 



X 
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r/2(l + V3), j/ii = -l + V3. (65) 



Using them to compute {x + yrf'), which is proportional to c, we get 

(x + yr?2)i0'^ = 277^, {x + yrj'f'''' = 

{x + yri^f^ = 27/2^3, (66) 

hence 

V3{X + 7/r?2)10a _ ^ y^2)ll ^ o_ (g7) 

Thus if hu = —s/^hioa, then x + yr) is zero, giving c = and the model considered in the last two 
sections. As discussed in the last section, m3 = follows automatically in this case. 
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2. Majorana Mass Matrix Mjv 

There is no difficulty in computing tfie Majorana mass matrix Mn as we shall see, but M^^ is 
needed for Mi,, so a matrix inversion has to be performed. Since Mjv and M^^ satisfy the same 
symmetry constraint, (24) and (28), their solutions must have the same form. Henceforth we will 
use the form of and the same letters with a prime to parametrize M^. In the case of (33), 
the relation between these two sets of parameters are 



Prom Table 1, we see that IR3xIR5=IR3 if B=l, 2, or 3. The residual symmetry Gb for S = 3 
is given in (32), and the other two, Gi = —G2 = +1, are given in Table 1. The alignments are 

{(f)^) = (1, 1)"^, {(j)^) = 1, and (0^) = 0, so only B =IR1 and IRS contribute. The Majorana mass 
matrix M^^ of the form given in (33), with P' = hi and Q' = —h^. The parameters P, Q for M^^ 
can be obtained from (68). 

B. p = 2, C=IR1+IR2 

1. Dirac mass matrix M^ 

From Table 1, we know that IR5xIRll contains IRl, and IR5xIR10 contains IR2, so once again 
B in (55) consists of IRIO and IRll. The corresponding representation Gb and the alignments 
are the same as those given in Sec. VA. The Dirac mass matrix Mj/ is the sum of and M^, 
obtained using respectively alignments {(f)^^"') and {<p^^'') for IR2, and alignment (^^^) for IRl. Mi, 
can be calculated from (55) and (62). It is of the form in (37), with 



Prom (39), we see that in order to reproduce the model with c = and / / 0, we need to have 
y = and hence hioa = 0. 

2. Majorana Mass Matrix Mjv 

Since IRlxIRl=IRl and IR2xIRl=IR2, if we let the two Yukawa coupling constants be hi and 
h'l respectively, then the Majorana Matrix M^^ is given by (37), with P = h^^ and Q = h'f^. 



P = P'/Di, Q = -Q'/Di, Di:=P'-Q'\ 



(68) 




(69) 



28 



C. p = 3, C=IR3+IR1 



1. Dirac mass matrix M^, 

As remarked in Sec. IVB3, the here is a composite, with the first two rows given bt the 
in VA, and the last row given by the first row of the in VB. The parameters x, y, z, w in (40) 
can now be copied from (64) and (69) to be 

z = -^hiQb, w = ---^/in?7^ (70) 

where x^'^"', x^^,y^^ are given in (65), hs is the Yukawa couphng to IRS, and h'^ is the Yukawa 
coupling to IRl. 

Again, as in VA, because of (67), to render c = we need to have hu = —V^hioa- 

2. Majorana Mass Matrix Mjv 

Both Mn and M^^ are of the form (40), with their parameters related by 

P = {P'S' - R'^)/Dsa, Q = -{Q'S' - R'^)/D^a, := (P' - Q')Ds,, 

R = -R'/Dsb, S={P' + Q')/Dsb, := {P' + Q')S' - 2R'\ (71) 

The parameters for Mjv can be partially copied from VA and VB: 

P' = hi, Q' = -h3, R' = hyV2, S' = h'i, (72) 

where the couplings are respectively for IR3xIRl^IR3, IR3xIR3^IR3, IR3xIR3^IRl, and 
IRlxIRl^IRl. 

D. p = 3, C=IR3+IR2 

1. Dirac mass matrix 

As remarked in Sec. IVB4, the here is a composite, with the first two rows given bt the 
in VA, and the last row given by the second row of the in VB. The parameters x, y, z, w in (40) 
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can now be copied from (64) and (69) to be 



2^3 



z = --^hiob, w = -^h'lirj^, (73) 

where x^^, j/^^ are given in (65), hs is the Yukawa couphng to IRS, and h'^ is the Yukawa 

couphng to IRl. 

Again, as in VA, because of (67), to render c = we need to have hn = —\/3hioa- 
2. Majorana Mass Matrix Mj^ 

Both Mjv and M^^ are of the form (43), with their parameters related by 

P = {P'S' - R'^)/Dia, Q = -iQ'S' + R'^)/D4a, D^a := {P' + Q')D^b, 

R = -R'/D4b, S = {P'- Q')/D^b, := {P' - Q')S' - 2R'\ (74) 

The parameters for Mn can be partially copied from VA and VB: 

P' = hi, Q' = -h3, R' = hy^/2, S' = h[, (75) 

where the couplings are respectively for IR3xIRl^IR3, IR3xIR3^IR3, IR3xIR3^IR2, and 
IR2xIRl^IR2. 

VI. NEUTRINO MIXING OF A (600) 

A parallel was mentioned in the Introduction between symmetry considerations of the old and 
the new mixing data. The zero reactor angle and the maximal atmospheric angle of the old data 
can be explained by a S3 symmetry, a group generated by the matrices F = diag(l, a;, u;^) and 
an order-2 unitary matrix G with an invariant eigenvector (0, 1, —1)'^. In order to have a chance 
to explain the solar angle by symmetry as well, we need another order-2 unitary matrix G' which 
commutes with G. Such a matrix with an invariant eigenvector (1,1,1)^ does exist and is contained 
in the group 54 D S3. These two invariant eigenvectors together give rise to tribimaximal mixing 
and the correct solar angle. 

With the new data of a non-zero reactor angle and possibly non-maximal atmospheric mixing, 
we saw in the last three sections that they can now be explained by the symmetry group A(150), 
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a group generated by F and G of (25) and (27). Again, to have a chance to explain the solar angle 
by symmetry, we need to find another order-2 unitary operator G' which commutes with G. The 
smallest group containing such an operator as well as the subgroup A(150) is A(600). 

This motivation for the full-symmetry group A (600) is based completely on symmetry consid- 
erations. There is another motivation based on dynamics which also has a parallel with the old 
data. The most popular dynamical model for the old data is based on the symmetry group 
[16] , generated by F and & above. If its neutrino coupling to the 1' and 1" valons are dynamically 
suppressed, then the solar angle comes out correct, the matrix G becomes also a symmetry, and 
the group is promoted to S^. 

The analogy of with the new data is A(150). By setting the parameter c dynamically to 
zero, we saw in the last three sections that the solar angle can also be explained. In the GAP basis 
where the generators of A(150) are the F, G of (25) and (27), the mixing vector for the model with 
c = is (0, 0, 1)^. It turns out this is just the invariant eigenvector of the operator G' in the group 
A(600). 

To see that, let us first find out what G' is. It must be an order-2 element of A(600) that 
commutes with G of (27). In terms of the generators of A(600) given in Sec. IIB, G = fic^cf. 
According to (7), /i, 03,04 all commute with the order-2 element ^3, so a natural candidate for 
G' is a three-dimensional representation of ^3. From Table 5, there arc 18 three-dimensional 
representations. We must choose one for / and g which reproduces F and G in (25). This turns 
out to be IR14. The representation of ^3 in IR14 is G' = diag(— 1, —1,1), which does possess the 
desired invariant eigenvector (0,0,1)-^. 

We know from Sec. IIB that F and G' generate a A4 subgroup of A (600), and we already know 
that F and G generate the subgroup A(150). In this way A(600) for the new data is like S4 for 
the old data, which contains and are generated by the matrices F, G, and G'. In fact, if we use V 
of (12) to convert the F-diagonal to the GAP representation, F is the same as F, G' is the same 
as G', though G and G are different. F, G generate an ^3 subgroup of ^4, and F, G generate an 
A(150) subgroup of A(600). 
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A. Left-Handed Mass Matrices for A(600) 



Apply the symmetry constraints 

F^M^F = M^, G^MyG = M^, G'^M^G' = M^ (76) 

to the Hermitian Mg and the symmetric M^,, we get the same solution Mg as shown in (26). For 
M,^, the extra symmetry constraint provided by G' reduces the form in (26) to the form in (18) for 
the special model. Hence the c = model of A(150) can be realized by the increased symmetry of 
A(600), just like the zero 1', 1" model of A4 can be realized by the increase symmetry of ^4. 

VII. SUMMARY 

The symmetry group A (150) is applied to neutrino mixing to obtain the correct reactor and 
atmospheric angles. Three methods are used to implement this symmetry: by imposing it to the 
left-handed fermions alone, by imposing it to both the left-handed and the right-handed fermions, 
and by constructing dynamical models based on the symmetry. The relative merits of these three 
methods are discussed. Generally speaking, as far as mixing parameters are concerned, it is suf- 
ficient to use the left-handed symmetry. A simple model with the parameter c = has been 
discussed in some detail. This model gives the correct solar angle and can be reproduced from the 
enlarged group A (600) by symmetry alone. 
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Appendix A: Real Representations of A(150) 



We see from Table 1 in Sec. II that the characters of IRl, 2, 3, 12, 13 are real. As noted in 
Sec. IIA, these are actually real rather than quaternionic representations because their Frobenius- 
Schur indices are +1. Let fi{n) denote the representation of the generator fi in IRn. For IRl and 
IR2, the irreducible representation is the same as the character, which is explicitly real. For IR3, 
the representation given in GAP 

c :)•'-(::.)• 



/i(3) 



/I 0\ /I 0~ 

/3(3) = f J, /4(3)=( ^ 

is not all real, and neither is the case for IR12 and IR13: 
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/3(13) 



/?7 0\ 

7?^ 

1 

1 

77"* 

Vo ??/ 



/4(13) = 



/?7 \ 

??2 



7]'^ 



r/'' 







(A3) 



\0 TiV 



Although these matrices may not be real, they are all unitary. In that case [11], there is 
a unitary and symmetric matrix B{n) which renders B{n)fi{n)B(n)~^ = fi{n) explicitly real. 
Moreover, S{n)fi{n)S{n)~^ = fi{n)* when S{n) = B{nf. Our task is to find S{n) and B{n), and 
here is how. 

To simplify writing, let us drop the dependence on n and i. If is an eigenvector of / with 
eigenvalue A, then v* is an eigenvector of / with eigenvalue A*. Prom BfB~^ = / = /*, we 
conclude that u := B~^v is an eigenvector of / with eigenvalue A, and w := B~^v* is an eigenvector 
of / with eigenvalue A*. If A is complex and non-degenerate, then v = Bu = (Bw)*, hence 



Su = w* 



(A4) 



because B is unitary and symmetric. Even if A is degenerate or real, (Su)* is still an eigenvector 
of / with eigenvalue A. Together with (A4), these relations can be used to obtain S{n) from Ui{n) 
and Wi{n). 

Once S is known, B = s/S can be computed as follows. Suppose is the normalized eigenvector 
of S with eigenvalue (Tq, then 

S' = ^ (Taeael, ^ S = ^ Voaeael- (A5) 

a a 

Since each has two values, there can be many iJ's, but any one of them would do. 



1. n = 3 

Only 72(3) in (Al) is complex, u = (1,0)-^ is its eigenvector with eigenvalue A = a;, and 
w = (0, 1)^ is its eigenvector with eigenvalue A* = w^. Thus Su = w* implies 

3(3) = r J =eie{-e2el (A6) 
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with ei = (1, 1)^/^2 and eg = (1, -1)^/^2. Hence 

(X X* \ 
X* X J 

where x = {l + i)/2 = exp{m/A)/V2. The real representations S(3)/i(3)S(3)-^ = /^(S) are given 

by 

/i(3) = , /2(3) = ' , 

73(3) =7M=(^ ), (A8) 



1 



where 



Cr := cos(27rir), := sin(27rir). (A9) 

We computed S{3) using (A4) to illustrate the general technique. In the present case, it can be 
obtained without this heavy machinery. Since /2(3)* is obtained from /2(3) by interchanging the 
(11) and the (22) entries, it is obvious that 5(3) must be given by (A6). 



2. n = 12 

Since the complex matrices /3(12) and /4(12) are diagonal, it is again easy to obtain S'(12) 
directly as follows. /3(12)* is obtained from /3(12) by interchanging the (11) entry either with the 
(22) or the (55) entry, but comparing /4(12)* with /4(12) tells us that it must be the (55) and not 
the (22) entry. From we then know that (22) must swap with (66) to get /|, then from we 
know that (33) must swap with (44), thereby obtaining 
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Taking the square root of -S'(12), we get 



S(12) = 
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fiil'i) arc then given by 
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a. n = 13 



The computation is similar to the n = 12 case. The results are 
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